The Flux-Across-Surfaces Theorem 
under conditions on the scattering 

state 

D. Diirr\ T. Moser\ P. Pickl^ 

1 Mathematisches Institut der Universitat Miinchen, 
Theresienstr. 39, 80333 Miinchen, Germany 

e-mail: duerr@mathematik.uni-muenchen.de, moser@mathematik.uni-muenchen.de 

2 Mathematisches Institut der Universitat Tiibingen, 
Auf der Morgenstelle 10, 72076 Tiibingen, Germany 
e-mail: peter.pickl@um-tuebingen.de 

February 5, 2008 



Abstract. The flux-across-surfaces theorem (FAST) describes the outgoing asymptotics of the quantum 
flux density of a scattering state. The FAST has been proven for potential scattering under conditions on 
the outgoing asymptote ^out (and of course under suitable conditions on the scattering potential). In this 
article we prove the FAST under conditions on the scattering state itself. In the proof we will establish 
also new mapping properties of the wave operators. 



1 Introduction 

The flux-across-surfaces theorem (FAST) is basic to the empirical content of scattering theory. 
The FAST describes the relation between the integrated quantum flux density of a scattering state 
over a (detector) surface and a (detection) time interval and the momentum distribution of the 
corresponding outgoing asymptote ■(/'out ■ In this paper we deal with the simplest case of one-particle 
potential scattering. We remark that the extension of the FAST to many-particle scattering theory 
is problematical, see |14) . 

With the quantum flux density (* denotes the complex conjugate) 

j"'' = Im(i/;*V^) 

and without spelling out the conditions under which it can be proven, the FAST reads 



\im J j j'^{x,t)-d(Tdt^ YnjY j j j"^' {x,t) ■ da dt ^ J \iiont{k)\'^d^k, (1) 

T RS T RY. Cs 

where S C 5*^ is a subset of the unit sphere, i?S] -.^ {x E M.^ : x = Rlj. a; e S} is the spherical 
surface covering the solid angle E and Cs := {fc e M'^ : ej, g E} is the cone given by S. Furthermore 
~ denotes the Fourier transform and V'out the outgoing asymptote to the corresponding scattering 
state -0 = ri+'f/'out with the wave operator $1+. 

The left hand side is interpreted and also shown to be the crossing probability of the particle 
crossing the surface i?S 0[71|2n|j [HllHl^]' From the crossing probability one derives the scattering 
cross section |lllll2j . The right hand side of Q relates the crossing probability to the S'-matrix. 
Technically, the FAST (dJ) has been proven requiring conditions on ^out • But it is clear that when 
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all is said and done one needs the conditions on the scattering state for which the FAST holds. 
In particular, the microscopic derivation of the cross section needs the FAST under conditions on 
the scattering state JT]. In the present paper we establish the FAST under conditions on the 
scattering state. 

The FAST has been put into a mathematically rigorous setting by Combes, Newton, Shtokham- 
mer in 1975 [7]. In 1996 the FAST was proven by Daumer et al. [S] for the Schrodinger case without 
a potential. One year later Amrcin, Pearson and Zuleta proved the FAST for short and long range 
potentials using methods in the context of Kato's i7-smoothness, requiring an energy cutoff on the 
outgoing asymptote 13^. (More precisely, suppV'out is compact.) This at first sight innocently 
looking requirement seems however to be an important hindrance towards the physically relevant 
formulation of the FAST with conditions on the scattering state itself. We shall discuss this in 
somewhat more detail later. In 1999 Teufel, Diirr and Berndl gave a proof based on eigenfunc- 
tion expansions without an energy cutoff, but still requiring smoothness properties of the outgoing 
asymptote for potentials falling off faster than x~'^ Panati and Teta gave a proof for the 

special case of point interactions under conditions on the scattering state j21) with similar methods 
as in "25~. In 2003 Nagao 19 proved a weaker result, namely leaving out the second equality in 
equation This proof works for short range potentials falling off faster than the dimension of 
the space {— 3) and requires only conditions on the scattering state. By leaving out the second 
equality in |^ the result does not establish the connection to empirical data of a typical scattering 
experiment, as it does not establish the probabilistic meaning of the quantum flux as a crossing 
probability or in technical terms it does not establish that the flux points asymptotically outwards. 
In the same year Diirr and Pickl '13' proved the FAST for a Dirac particle under conditions on the 
scattering state alone using eigenfunction expansions. 

We provide now a proof for the Schrodinger case combining the techniques of the proofs in 
|13l 1^ to establish the FAST under conditions on the scattering state and for potentials falling 
off faster than The idea is to prove the FAST under almost optimal conditions on 4'out, 

which can be translated to reasonable and easily checkable conditions on the scattering state. It is 
clearly essential that there is no energy cutoff on ■(/'out, because it is highly unclear whether there 
are any reasonable conditions on the scattering state ensuring a cutoff on ipout (cf. ©)■ Having 
formulated the task to prove the FAST under conditions on the outgoing asymptote which can 
be transferred to conditions on the scattering state we like to remark, that there are no suitable 
assertions in the literature which allow to transfer conditions on ip to ipout in the context of the 
proof of the FAST.^ We shall elaborate this further considering eigenfunction expansions. We 
recall the generalized Fourier transform (see Lemma which maps the scattering state ip to the 
ordinary Fourier transform ^out of ^out ■ 



where (p*^{x,k) are the generalized eigenf unctions. In Lemma El we collect the properties of the 
eigenfunctions which are — concerning smoothness and boundedness — in general very poor: The 
generalized eigenfunctions are solutions of the Lippmann-Schwinger equations: 



in which we note the appearance of the absolute value fc of fc in the spherical wave part. Derivatives 
of k of higher order than one behave singular for A; ^ 0. Therefore we expect in general that the 
derivatives of the generalized eigenfunctions (of higher order than one) are unbounded for small 
k.^ In view of (0) this singular behavior is typically inherited by i/iout and it is hard to see how 
"extreme" conditions on -^out like i/'out in Schwartz space or '0out compactly supported can be 
satisfied by reasonable scattering potentials or states. This caveat applies to the above cited works 
on the FAST except |E1 ^1 1^ • Our task is thus to read from (0) proper conditions on V'out 

^For mapping properties between ip and Vout, which are not applicable in our case, see e.g. 1271 . 
^For a point interaction the generalized eigenfunctions can bo explicitly computed 0, p. 39. They have exactly 
this singular behavior. 




(2) 




(3) 
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which can be formulated in terms of the scattering state and then to prove the FAST under these 
conditions. 

The paper is organized as foUows: In Section 13 we recall the mathematical basics of scattering 
theory including recent results and fix notations, in Section 13 we formulate and prove the FAST 
under weaker conditions on the asymptote than in |25| . The conditions will be transformed by the 
mapping Lemma|21to sufficient conditions on the scattering state. The most complete statement is 
Corollary^ Technically the FAST is proven by stationary phase methods, which turns out — due 
to our necessarily weak conditions — to be a rather involved modification of standard results, e.g. 
Theorem 7.7.5 in 15 . The proof of the modified assertion is done in the appendix. 

2 The mathematical framework of potential scattering 

We list those results of scattering theory (e.g. |3 UHl UEl dl EH 123 El ESI ) which are essential for 
the proof of the FAST in Section El 

We use the usual description of a nonrelativistic spinless system by the Hamiltonian H (we use 
natural units fi = m = 1): 

H ■.= -^A + V{x)=:Ho + V{x), 
with the real- valued potential V £ {V)n, defined as follows: 

Definition 1. V is in {V)n, n=2,3,4,---, if 

(l) L2(R3), 

(ii) V is locally Holder continuous except at a finite number of singularities, 
(Hi) there exist positive numbers e, Cq, i?o such that 

\V{x)\<Co{x)-''-' for\x\>Ro, 

where (•) := (1 + {-f)^ . 

Under these conditions (see e.g. ^HI) H is self-adjoint on the domain D(_ff) =D(_ffo) = {/ £ 
L^(IR3) : / |fc^/(fc)pd'^fc < cx)}, where / := J^f is the Fourier transform: 

f{k) := (27r)-i j e'^'^-^ f{x)d?x. (4) 

Let U{t) = e~'^*. Since H is self-adjoint on the domain D(i/), U{t) is a strongly continuous 
one-parameter unitary group on L'^iM?). Let (j) £D{H). Then (j)t = U{t)(f> gD(iJ) and satisfies the 
Schrodinger equation: 

d 

'''glMx) = H(f>t. 

We define the wave operators il± with the range Ran(il-i-) in the usual way: 

n± : L^{R^) Ran(r2±), 
n± :=s- lim e'"'e~""°\ 

where s-lim denotes the limit in the L^-sense. Ikebe [T^ proved that for a potential V G {V)2 the 
wave operators exist and have the range (this property is called asymptotic completeness) : 

Ran(r!±) =7ieont(i^) =Wa.c.(^r), 
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where HcontiH) and Ha.cXH) denote the subspaces of L^{R^) consisting of states, which belong to 
the continuous and the absolutely continuous part of the spectrum of H. Then we have for every 
"ip € Ha.c.i^) asymptotes ■0in,V'out £ ^^(M'^) with: 

fl-^pin ^ -ijj ^ ^ + i'out- (5) 

On T){Ho) the wave operators satisfy the so called intertwining property 
On 7^a.c.(i?)nD(i/) we have then that 

Hoflg^ = ng^H. (6) 

We will need the time evolution of a state ip g Tia.c. {H) with the Hamiltonian H. Its diagonalization 
on 7ia.c.{H) given by the eigcnfunctions Lp±: 

{-^A + V{x))ip±ix,k)^'^ip±{x,k). (7) 

Applying (-iA - ^ T M)^^ in one obtains the Lippmann-Schwinger equation. We recall the 
main parts of a result on this due to Ikebe in which is collected in the present form in j25| . 

Lemma 1. Let V E {V)2- Then for any k E M'^\{0} there are unique solutions (p±{-, k) : M.^ ^ C 
of the Lippmann-Schwinger equations 

Ztt J \x — x'\ 

with the boundary conditions \im^x^^^{ip±{x, k) — ^^i-k-x-^ _ yjhich are also classical solutions of 
the stationary Schrddinger equation Q), such that: 

(i) For any f G L'^{M!^) the generalized Fourier transforms^ 

(J^ifm = -L^ I i.m. [ ^l{x,k)f{x)d^x 

exist in L2(R3). 

(a) Ran(J^±) — L^(R'^) and J-± : T-La.c.{H) L'^{M?) are unitary and the inverse of J^± is given 
by 

{Tg'f){x) = -^\A.m. I ^±{x,k)f{k)d'k. 
(ztt) 2 J 

(Hi) For any f E L^(M'^) the relation fl±f = J-^^J-f hold, where T is the ordinary Fourier 
transform given by 

(iv) For any f E D(i7) HTia.c. (i?) we have: 

and therefore for any f E Tia.cXH) 

e'^"'fix) = [:F^'e~^'^'T±f) {x). 



'l.i. m. J denotes s- lim J . 
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In order to apply stationary phase methods we will need estimates on the derivatives of the 
generalized eigenf unctions: 

Lemma 2. Let the potential satisfy the condition {V)n for some n > 3. Then: 

(i) ip±{x, •) G C"^^(M'^ \ {0}) for all x £ M.-^ and the partial derivatives'^ d^(p±{x,k), 
\a\ < n — 2, are continuous with respect to x and k. 

If, in addition, zero is neither an eigenvalue nor a resonance of H , then 

(a) sup \ip±{x,k)\ < CO 

and for any a with \a\ < n — 2 there is a Ca < oo such that 

(iii) sup |Kl"l-ia^.^±(a;,fe)| < c„(x)l"l, withK:=^. 
fee R3\{o} ^ ' 

Similarly, for any /e{l,...,n — 2} there is a ci < oo such that 

(iv) sup -^ip±{x,k) <ci{xy. 
fee R3\{o} 

Remark 1. Zero is a resonance of H if there exists a solution / of Hf = such that {x)~^ f € 
L^(R^) for any 7 > ^ but not for 7 = 0.^ The appearance of a zero eigenvalue or resonance can be 
regarded as an exceptional event: For a Hamiltonian H — Hq + cV, c € K, this can only happen 
in a discrete subset of R, see pQ, p. 20 and 17 , p. 589. 

Remark 2. Lemma 12 except the assertion (iii) was proved in |25|. Theorem 3.1. Assertion (iii) 
repairs a false statement in Theorem 3.1 which did not include the necessary factor, which 

we have in (iii). For \a\ = 1 which was the important case in that paper there is however no 
difference. For completeness we comment on the proof of this corrected version in the appendix. 
We note that the problem which we address here comes from the appearance of the absolute value 
of k in the Lippmann-Schwinger equation |(HJ), see also the introduction. In fact Lcmma|21(iii) is to 
our knowledge the best one can say about the derivatives of the generalized eigenfunctions w.r.t. 
the coordinates. Note that the higher derivatives (|a| > 2) become unbounded for small fc. In 
it is claimed that the derivatives stay bounded for small fc, see Proposition 3.8. therein. The proof 
of this proposition apparently disregard the behavior of the coordinate derivatives of fc. 



3 The flux-across-surfaces theorem 

The FAST Q is a relation between a scattering state and its corresponding asymptote. As already 
emphasized it is important to establish the FAST with conditions only on the scattering state 
(and the potential V). Since tp — Q+'tpont we get by the well known expansion Lemma Q (ii)-(iv): 

iIj{x, t) = JF^^e~*"^*'(/'out(fc) and we can express the flux in by its asymptote '0out(fc)- Therefore 
we will proceed in the following way: We will first prove a FAST under conditions on i/^out {k) and 
then we will translate these conditions to the corresponding scattering state. The connection 
between a scattering state and its corresponding asymptote is given by the expansion Lemman(ii) 
and (iii), cf. ©.^ That means, as already emphasized in the introduction, that the properties of 
V'out are via (O inherited by the properties of the generalized eigenfunctions, which are in general 

*We use the usual multi-index notation: a = {01,02, as), Oj S No, |o| := oi + 02 + f*3 and d^f{k) := 
^There are various definitions, see e.g. 1271 . p. 552, 0, p. 20 and 1171 . p. 584. 

^Because of Lemma|^(ii) we can use the generalized Fourier transform without the 1. i. m. whenever £ 
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very poor, see Lemnia|21 especially (iii). More precisely, we will see later in the mapping Lemma 
13 that the decay properties (i.e. for large fc) of V'out(fc) and its derivatives depend mostly on the 
differentiability of ip{x), while the properties of V'out(fc) and its derivatives for small k are closely 
related to the corresponding properties of the generalized eigenfunctions ip*^_{x,k). Therefore we 
now define a class of asymptotes, , for which we can prove the FAST and which has the same 
poor properties for small k as the generalized eigenfunctions in Lemma [21 The exponents which 
determine the decay for large k are optimized to get a large class and are of technical interest. The 
class of the suitable asymptotes is defined as follows: (In the following definition we have the 
Fourier transform of V-'out = fl'^^^ (cf. jSJ) in mind.) 



Definition 2. A function / : \ {0} ^ C is in if there is a constant C G M+ with: 

l/(fc)l<c(fc)-l^ 

Kf{k)\<C{k)-^, \a\ = l, 



\Kd^f{k)\<C{k)-\ \a\=2, K = 



< c{ky 



With that class we can formulate a FAST under conditions on V-'out (fc): 

Theorem 1. Let the potential satisfy the condition (V)^ and let zero be neither a resonance nor an 
eigenvalue of H. Let tliout{k) S Q'^ ■ Then i{}{x,t) — e~*^*f2+V'o«t(a') is continuously differentiable 
except at the singularities ofV and for any measurable E C 5^ and any T £ M; 



lim / / {x,t) ■ dcrdt — lim 

R — >oo J J R — >oo 

T fiS T RY, 



r{x,t)-da- dt = j |Vo«t(fc)rrf'fc, (9) 



where := {x e : x = Ruj, a; G E} and Cs := {fe G : ej. G E}. 

The crucial condition in Theorem ^ is V'out(fc) G ■ We introduce now the class Q of scattering 
states for which we can prove that the corresponding asymptotes are in . 



Definition 3. / is in if'' 



G L2(R3), n G {0, 1, 2, 8}, 
G l2(M3), nG {0,1,2,3}. 



Then Q -.^{j e-'"'g°. 



That means G is a subset of 7ia.c.(-ff) and is invariant under finite time shifts, i.e. if / G ^ then 
^-iHtj ^ ^ Furthermore G is dense in 7ia.c.(-ff) which can be seen e.g. by the results used 

in 121, p. 5368: Let V4 := {g{H){x)-'^ilj\g G C^(]0, oo[), V' S L^{R^)}. Since our potentials have 
no positive eigenvalues (cf. Section EJ we have that T>4 C TLa.c.{H). It is easy to check that is 
dense in H^.dH). Moreover (cf. g]) we have that X»4 Q D{H) n D((a;)'^). Again by T" HVi C V4 
which implies that T>4 C Q. Hence, G is dense in 'Ha,.c.{H). Note that the condition ip G G can be 
formulated also more explicitly (cf. Remark 13). We wish to remark that the condition ip G C^{H) 
seems to be natural: Wave functions in thermal equilibrium are typically in C°°{H), see j26| . 

With Definition 13 we can state now the important mapping lemma. 



n=l 
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Lemma 3. Let V G {V)i and let zero be neither a resonance nor an eigenvalue of H . Then: 

i;{x) eg^ n^{k) = i^outik) e g+. 

The proof is adapted from ll^l and can be found in the appendix. The lemma holds also for f2_|- 
replaced hy fl- and V'out by ipin-^ 

Theorem n] and Lemma 13 give the following corollary, the FAST under conditions on the scat- 
tering state. 

Corollary 1. Let V € {V)4 and let zero be neither a resonance nor an eigenvalue of H. Let 
^ G G- Then for any measurable S C S"^ and any T e K; 

C30 C30 

^lim J J j'^{x,t)-d(Tdt^ Van j j j"^ {x , t) ■ da dt ^ J \$out{k)\'^d^k. 



T RT. T RT. Cs 



Remark 3. Instead of the condition ip £ Q one can give of course also the condition on ?/; and V 
more explicitly. In the following we will give two examples for -0 and V such that %Ij G . The set 
of wave functions Q for which the FAST holds is then — according to Definition O — given by the set 



teR 

Let H"^'^ the weighted Sobolev space 

Then one can find for example the following conditions for which ip £ Q'^ . 

(i) V e {V)2, V e C^'^{S? \ £), where £ denotes the set of singularities of V and 

(ii) V e {V)2, V e H^'^-^ n H^ * and i; e U^AH) n n ijS'l 
Clearly both sets for ip are dense in T-Ca.c.iH). 



Proof of TheoremQ] We will prove the flux-across-surfaces theorem (jSJ for some T > 0. This is 
sufficient since: (f < 0, T > 0) 



RS T RY. 



lim / / r{x,t) ■ dcrdt = lim / / ^ {x,t) ■ dcrdt, (10) 



with (in the second line we use Lemma (ii)-(iv), again without the 1. i. m., because of Lemma 12 

(ii) ami iioutik) e g+ C LH 



'iP{x,t) ='iJj{x,t + T -T) ^ {2Tr)-^ J e-^^e''^^^ipoutik)ip+{x,k)d^k 

=:(27r)-i j e~''^XoAk)^+{x,k)d^k. (11) 

It is easy to check that Xout(fe) € , if V'out(fe) G , which means that 5+ is invariant under 
finite time shifts. Hence, With pU|) and pi|l we get: 

oo oo 

^lirn^ y" J j^'ix,t)-d(Tdt^\im^J J j'^ (x.t) ■ dadt ^ J \xont{k)\^ d"" k ^ y |^out(fe)|'rf'fc. 



f RS T RT. Cs Cs 



*It would be interesting to have similar mapping properties for Q,^. 
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Of course, this argument is also valid for the integration over \j^{x, t) ■ dcrj 
Let T > be fixed. Using Lemma ^ (ii)-(iv) and we get: 



(2n) 



' ^out{k)(p+{x, k)d^k 



-^out(fc)e^'^-"d^fc + (2^)-^ / e-'—^Pontik)rjix,k)d-'k 



= : a{x, t) + (3{x, t). 
The flux generated by this wave function is: 

/•^(a;, t) = Im(a*Va + a*V/3 + (3*Va + ^*V/5), 



(12) 



(13) 



where a is obviously continuously differentiable and for the differentiability of j3 see (20) and 
(28)-(30). In |5] and [53] the function a{x^t) is estimated using the formula 

a{x, t) = {27Ttt) I Vout(y)d'y (14) 

and conditions on '0out(a')- According to Lemma|31we can control ipontik), but not ^Jo^tix). Hence, 
we have to estimate a{x,t) directly in terms of ipoutik). This will be done by using stationary 
phase methods. First, we will calculate jg = Im(a*Va) by using Lemma^J which is formulated 
for a special class of wave functions IC D . This set has similar weak conditions as the set 
due to the necessarily poor properties of 'tpout{k) (see the discussion before Definition |2J|. Again we 
give here optimized decay properties, which are, however, not that strong as in the case of . 

Definition 4. A function / : R"^ \ {0} C is in K, if there is a constant C E M+ with: 

l,/(fc)l<c(fc)~^ 

\d^f{k)\ < C, \a\ = 1, \Kd^fik)\ < C{k)-\ |a| ^ 2, 



With that class of wave functions we can formulate 



92 



Lemma 4. Let x(fc) he in IC. Then there exists a constant L G M-|_ so that for all x G 



X{k)d'k~ ( ^ 



'X{ks 



< 



L 

t2' 



(15) 



where ks — 



The proof of the lemma can be found in the appendix. 

Applying that lemma on a (a;, t) in ^ we get, with an appropriately chosen constant L: 



a[x,t) - ( - 



^^^out (f ) 



< 



L 

i2 



and analogously: 



Va(a;, t) — i \ — 



(f ) '^-^ (f) 



L 



which gives for the flux j'g = Im(a*VQ;): 



< 



t2 



(16) 



(17) 



(18) 
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We begin with the first term in 113|l for times t > Re; (We choose R big enough, so that 
iit > T.) 



(19) 



Inserting the asymptotic expression p8|l for the flux Jq we get instead of H19|l : 



4 S 



i4 



dndt 



^out (fc) 



k^dfldk. 



(20) 



E 



where we substituted fe := (j^Oll gives in the hmit aheady the right resuh: 



Re 



Um 



S 



|V'out(fc)rA:^dMfc = J |i^out(fc)rd'fc. 



(21) 



From p8|) - (|20|) it is clear that also the modulus of Jq gives the right result. Hence, by justifying 
the use of the asymptotic expression for j'q , showing that the integral p9|l is negligible for times 
smaller than Rb (and large R) and by proving the smallness of the contributions of the three other 
terms in (|13|) we get Theorem ^ 

Using H18fl we can estimate the error between p9|l and (|20|) : 



L I I Rh~idndt^—R-^. 



(22) 



rI ^ 



which tends to zero for large R. 

We evaluate now the flux integral for times smaller than Re : 



R'S 



J J j^{nR,t) -nR^dfldt. 



(23) 



T S 



Substituting t ^ Rt we get: 

1 

R^'6 



R~s 



J J jt{Rn,tR) -nR^dfldt < J j \a{Rn,tR)\\V ^a{x,tR)\^=Rr,R^dVtdt. (24) 



z s 



We estimate a and Va separately. We start with a: 

aiRn, tR) ^ {2Tr)-i f e-''^^^^''^^4^out{k)d'^k 



(25) 



The exponent of the e-function has the stationary point at fcgtat = 7- Since t E [-^,R^i], fcgtat S 
, ^ [. Big momenta should be negligible, hence we divide the integration over k in small momenta 
up to fc < R^ and larger ones. This will be done by the following functions: 



1 for /c < ^i?^, 
/i(fe) = <j cos2 ((fc - ii?^) f ) for < /c < + 1, 
otherwise, 



(26) 
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'O for k < ^R^, 



/2(fc) = < sin^ ^ - f foi" <k< + 1 



(27) 



1 otherwise. 

We have then /i(fe) + /2(fe) ee 1 and get for if^ : 



,{Rn,tR) ^(271)--^ J e-'*(^^-'=T^)^^out(fc)/i(fc)d^fc 



^out(fc)/2(fc)d'fc=:/i+/2. 



(28) 



We choose now R large enough (such that ii?8 > 1), which means that the first integral in 
has no stationary point anymore. We will do two integration by parts: 



h - (2^)-i j e-^*('#^-'=^)V^out(fc)/i(fe)d'fc 

= (2.)-i I (v.e-*(^«-^^)) . r*(^^'^-^") 
= -(27r)-i /e-^*('#«-'='^) f Vfc • 
= : (27r)- 



-i{Rtk - Rn) 



Rtk-Rn\'^ 
e-^*('^^-'='^)(Vfc-g(fc))d3A: 



^out(fc)/i(fc)rf'fc 
^out(fc)/i(fc))) 



= (2^)-i y e-"(Vfl-fe^) ( Vfc 

The gradient can be written as: 
-i{R.tk - Rn) 



A{Rtk - Rn) 
\Rtk - Rn\^ 



(Vfc-g(fc)) ]d'k 



\Rtk - i?n|2 



(Vfc-g(fc)) - ^ 9fe 



-i{Rtkj — Rnj) 
\Rtk - Rn\^ 



(29) 



(30) 



A straightforward calculation yields for the right hand side of (|30|l : (we consider one summand) 



i?"^^|^out(fc)||/i(fc)| , ^ Rt\dkAi^out{k)fi{k))\ 

rrrr; ^^"t;; ^2 



\Rtk - i?n|4 |i?tfc - i?n|3 

i?i|afc^(^^out(fe)/i(fc))| ^ |9fc,9fc,(^out(fc)/i(fc))| 



|i?tfc - Rn\^ 



C4- 



\Rtk - Rn\ 



with constants Ck > 0, k — I, 2, 3, 4. Since < fc < ^Re + 1 and < t < R e we have: 



|i?ife - i?n| > -R- Ri > -R, 
I 1-2 - 3 



(31) 



(32) 



if R is large enough. Using and the definition of fi{k) we find, with an appropriately chosen 
constant M > 0, instead of H31|l : 



dk 



-i{Rtkj — Rnj) 
\Rtk - RnP 



(dkMk)) 



<-^|^out(fc)| 



+ ^ (l^^out(fc)l + \dk,Mk)\ + \dk^$out{k)\ 

V'out(fe)| + \dk,tpout{k)\ + |9fc;9fc^Vout(fe)| 
^ {\knt{k)\\dk^dkji{k)\) . (33) 



Mt 

M 

M 



3 THE FLUX-ACROSS-SURFACES THEOREM 11 



Using |9^^/)out(fc)| < C{k) ^, \a\ < 1, we get by and an appropriate constant M' > 
with: 

+ ^^^l^ j {k)-'\dk^dkji{k)\k^dkdn. (34) 

To integrate the second derivatives we use that {nd^'tp out < C{k)~'^, \a\ = 2 and 
k\dkidkj fi{k)\ < C{k). Hence, with an appropriately chosen constant C" we arrive at: 

We estimate now l2- Since V'out(fc) G we have: 

I/2I < {2t:)-^C J {k)-^''d^k < C"R^^, (36) 

with an appropriately chosen constant C" > 0. Hence, we find: 

\a{Rn,tR)\ =|/i +/2I < {C + C"){1 + tfR-^ =: C'{l + t)^R-^. (37) 

In a similar way we can estimate Va by: 

\^Mx,tR)URr.<C'{l + t)R-\ (38) 

To get this estimate we split again the analogous integral to H25|l into small and big momenta. The 
first part will be estimated by one partial integration using that \d^tpout{k)\ < C{k)~^, \a\ < 1, 
and lud^^poutik)] < C{k)~^, \a\ ~ 2, the second part (which is analogous to by using that 

IV'out(fc)| < C(fc)^i°. Inserting ^ and (gHl into ^ we get: 

1 1 

Re Re 

J J \a{Rn,tR)\\V^a{x,tR)\^^RnR^dndt < inC'^ J {l + t)^dt, (39) 

Z S 

R 

which tends to zero for R ^ 00. 

It remains to show that the three other terms in (|13() are neglig ible. In (Equations (15) 
and (16)) the function P{x,t) is estimated for some i?o > by: 

sup \Pix,t)\ < c ^ Vi? > 0, (40) 

sup |V/3(a;,i)| < c \ Vi? > i?o, (41) 
xe^H n[t + K) 

for t>T. The constant c depends on T, V'out(fc) and ^V'out(fc), and is finite for V'out(fc) G (cf- 
(20)-(28) in It is also shown that the last term in (O is neghgible (cf. p. 10 in [23). In ^ 

there are also estimates on the a(x,t) terms, but not under the conditions which we must require. 
We start with the second term in (fT^ : 

ao 00 00 

J J Im{a*V P)R^ndndt < J J \a\ |V/3| R^dQdt < J J \a\ -^j^^—^R^dndt. (42) 

T S T S OS 
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We divide again the time integration into two parts: 



E 



R{t + R) 



R^dQdt 



a 



s 



R(t + R) 



R^dVLdt 



R 



R{t + R) 



R^dVldt. (43) 



Hence, with (|37|l the first part is: 

5 

Re 

\a{Rn,t)\ 



R'e 



E 



R(t + R) 



R^dVLdt = / / \a{Rn,tR)\ 



s 
1 



< 



C'c{l + 1) 
R 



i?2(l+0 



dndt. 



R^dVLdt 



(44) 



E 



which tends to zero for R oo. 

It remains the second term in (|43ll . Applying the asymptotic expression H16|) for a we get: 



5 E 



l«(^»,^)l ji{t + R) '^^'^^- 



Ri ^ 



i?(t + R) 



dVtdt 



R.'B 



<- 



. Aire 

Vr 



i E 
1 

Re 



L cR^ 

R{t + R) 



dndt 



ipout (fc) 



1 _5 

^dfc + AttcLR 6 , 



(45) 



where we substituted k :— Since V^out G 5+ the bound in is finite and tends to zero for 
R ^ oo. The third term in can be treated anafogously to (|42|l - 145f) . ■ 



4 Appendix 



Proof of Lemma 12 Lemma |21 is proven — following the idea of Ikebe ^BJ — in [21]' The latter how- 
ever contains a mistake concerning the assertion (iii) , which overlooked the need for the smoothing 
factor K — , which puts the higher derivatives of the generalized eigenfunctions into the "right" 
Banach space. The need for this smoothing factor arises from the derivative of k which appears in 
the spherical wave part in ||SJ), see also the remarks in the introduction. Observing that, the proof 
goes through verbatim. Our statement (iv) follows also from the proof in |25| . replacing coordinate 
derivatives by the derivatives after k. In this case we note that there is no need for any smoothing 
factor. ■ 

Proof of Lemma Let iJj G Q. Then there is a x G and a < G R with: 

Using the intertwining property @ we get: 

V-out = = ^+'e-'"'x - e-'""'n^\ = e-^^"*Xout. (46) 
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Since is invariant under multiplication by e ' 2 * it suffices to show that Xout(fc) is in Q'^ . Let 
X e Since (a;)2iJ"x(a;) G i^(lR^), < n < 8 and (a;)''i?"x(a;) G L^{R^), < n < 3 we have: 



iJ"x(a;) G Li{R^) n L2(r3)^ < n < 8, 
(xyH^'xix) e ii(M^) n L2(r3)^ < n < 3, j = {1, 2}. 
Using the intertwining property © and Lemmas (ii), (iii) (cf. l(2Jl and Footnote we have: 

yXout(fc) = H^Xoutik) = T{H^ni\){k) = nVLl'Hxm 



(47) 



(27r)-4 / ^\{x,k){Hx){x)d^x. 



Applying n times on Xout(fc) (0 < n < 8) we get: 



^Xout(fe) = (2^)-t y ^\{x,k){H-x){x)d'x. 

Since the generalized eigenfunctions are bounded (Lemma El (ii)) and i?"x G Li 
we have with an appropriate constant C: 



IXout(fc)| <C(fc)-i^ <C7(fc) 



\-15 



(48) 

(49) 
< n < 8 

(50) 



Because of Lemma El (iii) and (I47|l we can differentiate Xout(fc) w.r.t. the coordinates and get an 
appropriate constant C with: 



\dkSont{k)\ 



(2^)-i / {dk,^X{x,k))x{x)d?x 



< C, Vfe e \ {0}. 



(51) 



Applying iJp three times in ()48|l and differentiating w.r.t. ki we get similarly to (|51|) : 
/c6afc,Xout(fe) = 8(2^)-i / (afc^(^;((r,fc))(i/3^)(a;)d3a;-6fc5xout(fc)^. 



(52) 

Again the right hand side is bounded because of Lemma El (iii), (I47|l and (|50|l . Hence, we get 
together with (|51|l : 

|afc.Xout(fe)| < Vfc e \ {0}. (53) 

To control a second derivative with respect to the coordinates we have to multiply by the factor k, 
since then the derivatives of the generalized eigcnfmictions Lp± are bounded by c(x)^, see Lemma 
El (iii). Hence by 



8(2^)-i {ndk,dk,^*+{x,k))x{x)d'x 



<C,ykeW^\ {0}. (54) 



\i^dk^dkSoMt{k) 
Using (|^ we get: 

k'^Kdk^dkSo.tik) =8(27r)-i j {ndk^dk^ifXix, k)) {H^x){x)d^x - 30k^'^jKXout{k) 
6k^'-^KdkjXout{k) - 6k^Xout{k)K 



^ Xout(fc) - 6fc^Xout(fc)K ^'^'^^, ^'^^ - 6fc^^K5fc,Xout(fc), (55) 

k k-^ k 



where the right hand side is bounded because of LemmaEl (iii), 1471) . H49|l and (|53|l . Hence: 

l«a^£out(fc)| < C(fc)-6 < C{k)-\ \a\ = 2,ykeM.^\ {0}. (56) 

implies also: 

\dkXout{k)\ < C{k)-^, Vfc e r3 \ {0}. (57) 

Applying Hq two times in ()48|l and differentiating two times w.r.t. k we get by LemmaEl(iv), (|47|) . 
(|5n|) and ||S7|) analogously to (|^ : 



|aIxout(fc)| < C{k)-^ < C{k)-\ Vfc e \ {0}, 



(58) 
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which means that Xout(fc) G G^- B 



Proof of LemmaEI At first sight Lemma 01 looks hke a standard stationary phase result, e.g. 
Theorem 7.7.5 in jl5j . But in our case we have (by necessity) very weak conditions on the function 
x{k), since we need to use the lemma for x(A;) = ?/'out(fe)- Especially the second derivative of 
x{k) w.r.t. the coordinates becomes unbounded for fc ^ 0. Furthermore, the stationary point kg 
is moving with x and t. 

First, we extract the leading order term of the integral (|15|l 



{X{k)~x{ks) + X{ks))d'k 



X{ks)d^k 



- {X{k)-x{ks))d'k. 



The leading order term can be easily calculated: 



X{ks)d'k = ( ^ 



(59) 



(60) 



We will now calculate the error between the left hand side of H59I) and the leading order term H60() : 



- {X{k)-x{ks))d-'k. 

The following splitting of the integration area turns out to be convenient (cf. Figure 1): 
Ai :={fc e : fc = |fc - fc,| < — }, := {fc G : fc < 2ks A |fc - k^] > — }, 



(61) 



^3 :={fc e M3 : fc > -fcj. 

The areas A2 and A3 have a small overlap. This is due to 
the use of suitable moUifiers. In Ai and A2 we shall perform 
two partial integrations w.r.t. the coordinates, in A3 we shall 
perform the derivatives w.r.t. k. Our proof will assume a; ^ 0. 
The case x = must be handled separately, but is much 
easier than the proof we give. It can be done by two partial 
integrations w.r.t. k similarly to our procedure which handles 
the area A3 (|^ . 

We first divide the integration area into Ai U A2 and A3 
by using the moUifier p(fe) : 



p{k) 



1, for k < |fcs, 



e exp 



V 



for |fcs < < 2fcs 



(63) 



(62) 




Figure 1: Sketch of the three inte- 
gration areas in the fe-frame. 



^0, for k > 2ks. 
The mollifier has the following properties: 

suppp = Ai U A2, \p{k)\ < 1, |1 - p{k)\ < 1, 



M 



There is an Af > such that: \dkp{k)\, \d^p{k)\ < — , |a| = 1, 
and\dlpik)\,\d^pik)\<^, |«| = 2. 



(64) 
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Using p we can write for l|61(l : 



+ 7' 

=:/i2 + Is- 

We start with the estimation of /12. We define: 

f{k) :=p(fe)(x(fc)-x(fcs)), k :^k k 
and get with two partial integration w.r.t. to A; : 



(i-p(fc))(x(fc)-x(fcs))rf'fc 



(65) 



(66) 



12 



1 

't 
1 

't 
1 

'"^ 
1 

1^ 



(Vfc 



^2 



/(fc)d3fc 



/ k-Vkf{k)~f{k) 



(Vfe£ 



/fc2 



ai| + |a2|=2 



<- 



/(fc)-fc-Vfc/(fc) 



fc4 



i2 



^ 5: r^r^9-9-/(fe) 



|ai| + |a2|=2 



(67) 



Because of the definition of p, the integration area in H67I) is Ai U ^2 (cf. H62|l and H63|) ). We will 
divide this area into Ai and ^2- Hence, /12 is estimated by: 



1^2 I < 



t2 



f{k)-k-VkI{k) 



Ai 
1 



fc4 



i2 



f{k)-k-Vkf{k) 



1 

" i2 



^ 5: r^r^arar/(fc) 



fc4 



ai| + l"2|=2 



d'^A; 



/c4 



|ai| + |a2|=2 



d'^fc 



A2 



f{k) - k ■ Vfc/(fc) 



k^ 



d^k + ^ 



A2 



^ 5: r^r^arar/(fc) 

^ |Qil + |a2l=2 



= -h+h- 



(68) 



We estimate /i first. With (|62|) and H63fl we see that for fe e Ai, p(fe) = 1 and thus we have for 
(|66|) : /(fc) = xC'') ~ x(fcs)- Using Taylors formula and then substituting fc by fc (cf. (|66|l 'l we get 
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for the first term ll of Ii : 



Ai 



fik) ~ k ■ Vkfik) 



fc4 



Ai 



x{k)-x{ks)-ik~ks)-VkX{k) 



1 

'1^ 



1 



Q!i| + I"2|=2 



2(fc-fes)4 

ai| + l"2|=2 



d-^k 



2k^ 



d^k, 



(69) 



where ^ is a vector between kg and fe. Hence we have ^ > -y- Using Definitional i.e. that 
dk.dkjXik) < Ck-^, we get for 



1 9C /■ 1 ,o~ 
^ ^ ' d^k< 



367rC 



ISttC 



(70) 



The second term of Ji can be estimated analogously: Instead of ^ we have k — k + kg with k > 
It follows that Ii is of order uniform in kg. The estimation of I2 is very similar, but p{k) 7^ 1 on 
A2. We use the volume factor d'^k integrated over ^2- Hence, it suffices to show that the integrands 
of the two terms of I2 are bounded by p- or by -p^ with some constant L > uniform in kg. The 
first integrand is: 



/(fc) - fc • Vfc/(fc) 



k^ 



< 



< 



p(fc)(x(fe)-x(fes)) 


+ 


|Vfc/(fc)| 


k^ 




/c3 



X{k) - x{kg) 



k^ 



E 



fc3 



E 



l(x(fc)-x(fcs))afc.p(fc)| 



fc3 



By mean value theorem there exists a ^ € (fcg, fe) with: 

Ix(fc) - x{ks)\ = |Vfex(OI \k - < CA: + Cfc„ 



(71) 



(72) 



since x(fc) G ^- Using (O, ((Sill, fc e ^2 (which means: k < 2kg, k > ^) as well as \dkiXik)\ < 
C, i = {1, 2, 3} we obtain: 



/(fe) - k ■ Vfc/(fe) 



fc4 



<- 



32C 16C 24C 48CM 24CAf 



jU3 ' 1,3 



fc3 



fc3 



(73) 



Similarly we estimate the integrand of the second term of I2 (I68|l . We pick one summand: (|q:i| + 
|a2| = 2) 



lk"'k"'d^^d^^f{k) 



< 



dl'dl-f{k) 



< 



4ix(fc)-x(fcs)ii5rgrp(fc)i , ^K'pmK'x{k)\ 

1.2 1.2 

s s 

, m'pik)\K'xik)\ , m'd^'xik)\ 



fc2 

8CM 4CAf 8CM 4C 
<^-^ + — ^ 4- + 



fc2 



fc3 



fc3 



fc3 2 



(74) 
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It remains to estimate I3 (|65() . We introduce a convergence factor Pe{k) 



, 1, for k < -, 



with < e < min(2p; 1). Then we get for I3 (gSJ): 

Is= [ (1 - p{k)) ixik) - x(fcs)) d'k 

e-^'T'+^^-^ (1 - p{k)) (1 - Pe(fc) + P.(fe)) (x(fc) - X(fcs)) d'fc 



= Ihn y (1 - p(fc)) (1 - p,(fc) + p,(fe)) (x(fc) - x(fes)) d^A: 

= lim 1 e-**+»'^- (1 - p(fc)) p,(fc) (x(fc) - x(fcs)) d'fc 

+ Ihn 1 e-^^t+^f'- (1 _ (1 _ p,(fc)) - x(fcs)) d^fc 

= Hm / e-*^*+^'^- (1 - p(fc)) p,(fc) (x(fe) - x(fes)) rf'A: 

+ Ihn 1 e-^^*+^'=- (1 - p,(fc)) (x(fc) - x{ks))d'k, (76) 
fc=i 

since 1 — p = 1 on supp(l — p^) (of. and ((TSJ)). The last term in the last Hne of (|7^ is zero 
(since x{k) G K, and by a standard Riemann-Lebesgue argument) and we get for 73: 

/3=lhny e-'*('^-'=-M {I- p(k))p,{k) {x{k)-x{ks))d'k 

=:lim / e~''i^~'"''') Mk,k^)k^dkdn. (77) 

We will perform now two partial integrations w.r.t. k: 



1 f fk^fe{k,k 



lim - / e '"v ' dkdn 

e-^ot J \k-ek-ksj 

j..^/.-(*-".)ft(^a.(^||*« 



<^ hm 



1 n f k^ feik, kg 

-Ok 



k Gk ' ks \ k Gk ■ k 



dkdQ 



hni J \D\dkdn 

<^hm J \D\dkdn + ^\im J \D\dkdn+^\im J \D\dkdn 



2k^<k<\ fe>i 

= :/i+/|+/|. (78) 
We start with the estimation of . Because of the integration area it suffices to show that D is of 
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order ^. Since Pe{k) = 1 in this area, D is given by: ((•)' denotes the derivative w.r.t. k) 



\D\< 



{k- Bk- ks)"^ 



|((l-p(fe))(x(fc)-x(fcs)))"| 



+ 



1 



1 



1 



|((l-p(fe))(x(fc)-x(fcs)))'| 

fc2 



.|(i-p(fe))(x(fc)-x(fe.))h 

We shall use (sometimes in a slightly modified version): 

fc2 fc2 



_ (fc fcg + fes)^ < q f I. > ^t. 
(fc - efc • fcs)2 - (fc - A:.)2 " (fc-fc.)2 1°'-'^- 2''^ 



< 



(79) 



(80) 



Using (IHOll we get instead of ifT^ : 



|/?|<9|((l-p(fc))(x(fc)-x(M))" 
47 



39 



(fc - fcs) 



2 I(l-P(fc))(x(fc)-X(fcs))|. 



|((i-p(fc))(x(fc)-x(fes)))'| 



(81) 



Using x(fc) € A^, i-e. 

|(x(fc)-x(fc.))'l<c(fc)"i<c, 
|(x(fc) - x(fc.))"| < c(fc)-2 <C{k)-' < C(l + fc.)-\ 

since fc > fcg, lEH), and (jHO)) we find that: 

818CM 9C 



1^1 <- 



1 + fcs 



(82) 



(83) 



It follows that /| is of order t ^ uniform in fcg . It remains to estimate /| and /| . First we consider 
"large" fcs : Let 2fcs > 1. Z? on the integration area of /| (where ^ > fc > 2fcs) is bounded by: (we 
use again |x'(fc)| < C{k)-^ and \x"{k)\ < C{k)-^) 



\D\< 



{k- Bk- ksf 



1 



1 



1 



\x(k)\ 



<- 



4C 



[k- Bk- fcs)' 

e V 

k ■— • fcg J k — Gfc ■ ks \k Gfc ■ ks J \k Gfc ■ k^ 
10c 52C 



Ix(fe) - x(fcs 



(84) 



^(fc)2 (fc-fcs)(fc) (fc-fcs)2' 

where we used H8()(l (and analogous estimates) with fc > 2fcs. Inserting H84(l into 177|l we get: 

\ll\^L\ini J \D\dkdn < ^ Win J \D\dkdn ^ ^ J \D\dkdn, (85) 

2fc,<fe<i 2fe3<fc 2k,<k 

which is integrable imiformly in fcs for 2fcs > 1. Hence, /| is of order uniformly in fcg for 2fcs > 1. 
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Similar we can estimate /| since then we have for D: 



\D\ <A\{p,{k){x{k)-x{ks)))" 
26 



10 



+ 



p,(fc)(x(fc)-x(fcs))| 



\{p,{k){x{k)-x{ks)))'\ 



<4 |(p,(fc)(x(fc) - x{ks)))"\ + 20 |(p,(fc)(x(fc) - x(fcs)))'| 
+ 104|p,(fc)(x(fc)-x(fcs))|. 



(86) 



The integration of H86|l over the area k > ^ yields a uniform bound in kg and e for 2ks > 1. It 
remains to estimate /| and /g for 2fcs < 1. /f can be estimated analogous to H86() since A: > i > 1 
and we have again: 



' <^<2. 



For /| we split the integration into: 



/g <— lim 



\D\dkdn + — lim 



(87) 



\D\dkdn. 



2fe3<fc<l 



l<fc<i 



The second term of the right side of H88|l can be estimated analogous to /| for 2ks > 1. Thus 
remains the following integral: 



7^ lim 

2fc3<fc<l 

The integrand is bounded by: (we use again (|Sn|) ') 



\D\dkdn. 



(89) 



\D\ < 



+ 



{k Gfc • fcg)^ 

fc2 



lx"(fc)| + 



1 



<4C 



D^l + \D2\ + iDsl. 



(x(fe) - x(fe.)) 

fc2 



lx'(fc)l 



■(x(fe)-x(fes)) 



(90) 



We have to integrate D over a bounded interval. Hence, Di yields a uniform constant. The 
derivative in D2 has at most two zeros in ^3. So we can divide the integration area into three 

subsets on which dk ((Fi^xi^) '^^^^ change the sign. Then we can apply the fundamental 
theorem of calculus to conclude that also the second term yields a uniform constant, using H80|) . 
Z?3 can be written as: 



■(X(fe) - x(0) - fcx'(O) + X(0) - x(fes) + kx'm 



■(x(fc) - x(0) - fcx'(O) + fcsg(fcs) + fcx'(O)) 



< 



k ' / k ' ^1 



■ (x(fe)-x(o)-fcx'(o)) 



+ 



--: \Dk\ + 



fc2 



1 



k ' / k ^k ' 

2 1 
3 ' 



■(fcs5(fes) + fcx'(0)) 



(91) 
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with appropriate bounded g{ks) : By Taylors formula and since |Vfex(fc)| < 3C we get 

\g{ks)\<3C. 

J5| in (|91|l can be treated analogous to D2 since the derivative has at most five zeros and 



<40C. 



(92) 



(93) 



To get l)93|l we use again estimates like (|80f) with k > 2fcs. Now we estimate in H91|) . Since the 
integration area is bounded it suffices to show that is uniformly bounded: 



1^3! < 



^— )'(x(fc)-x(0)-fcx'(0)) 



+ 



(94) 



Using Taylors formula we can linearize the x(fc)-terms and get: (0 < C < 1) 



Using |x"(fc)| < C and (|80|l (again we use also similar estimates with k > 2ks) one gets: 

\dI\< 120c. 

Hence, is of order uniform in fcg. It follows Lemma 0] ■ 



(95) 



(96) 
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